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We provide a simple transfer function that determines the effect of an early matter-dominated era 
on the gravitational wave background and show that a large class of compactifications of superstring 
theory might be tested by observations of the gravitational wave background from inflation. For 
large enough reheating temperatures > 10 9 GeV the test applies to all models containing at least 
one scalar with mass < 10 12 GeV that acquires a large initial oscillation amplitude after inflation 
and has only gravitational interaction strength, i.e., a field with the typical properties of a modulus. 



I. INTRODUCTION 

Our description of the forces of Nature consisting of 
the standard model of particle physics and general rela- 
tivity as the theory of gravity is astonishingly successful. 
However, the theory is incomplete. On large scales grav- 
itational dynamics requires the introduction of 'invisible' 
components like dark matter and dark energy. The stan- 
dard model of particle physics is lacking neutrino masses 
and contains many unexplained parameters. The most 
serious problem, however, is that we have no quantum 
theory of gravity. Classical general relativity has singu- 
larities. It breaks down when the curvature becomes too 
large. The most developed attempt for a quantum theory 
of gravity is string theory which may turn out to be the 
fundamental theory of Nature. 

String theory has the weakness that it is very difficult 
to test experimentally. This is not only due to the fact 
that 'stringy' effects become relevant only at very high 
energy, but is also a consequence of the landscape [T] . It 
has turned out that string models are extremely versatile 
and able to predict more or less everything. However, if 
we want to take some version of string theory seriously as 
a physical theory of Nature, it needs to be falsifiablc. In 
this paper we propose a test for the existence of at least 
one scalar field in the model that affects the cosmic evolu- 
tion after inflation. The well-known moduli problem [2] 
predicts a phase of matter dominated expansion some 
time after inflation. This matter has either to be diluted 
by a subsequent phase of 'thermal inflation' [3] or, more 
naturally, the moduli have to decay. In this paper we 
concentrate on this second possibility which does not re- 
quire any additional ingredients. Qualitatively, however, 
our discussion also holds for thermal inflation [4]. 

Moduli fields describe the configuration of the curled 
up extra dimensions of the string compactification. They 
have gravitational coupling strength only which makes 
them long-lived. They must be stabilized in order for 
the measured masses and couplings of standard model 
particles to have well-defined constant values as ob- 
served. Therefore, they must have finite mass. In su- 
perstring models their mass is typically of the order 



of the gravitino mass m 3 / 2) so that their decay width is 
r ~ m^/mpj ~ TOg^ 2 /mpj, where m p i denotes the Planck 
mass. If they are displaced from the origin after infla- 
tion, they perform oscillations and the Universe becomes 
matter dominated soon after. 

If the moduli were cosmologically stable they would 
overdose the Universe and if they would decay after 
^bbn ~ 0.1 s the success of big bang nucleosynthesis 
(BBN) were spoiled, because neutrinos would not have 
had enough time to thermalize [5]. Altogether, the mod- 
uli must decay before BBN. 

We show that such a matter dominated phase before 
BBN leaves a detectable imprint on the gravitational 
wave background from inflation: it significantly reduces 
the amplitude of the gravitational wave background at 
frequencies accessible to ground- or space-based detec- 
tors, compared to those probed by observations of the 
cosmic microwave background (CMB). Such a spectrum 
could be ruled out if a gravitational wave background 
from inflation were detected not only by future CMB ex- 
periments, but also by gravitational wave experiments at 
higher frequencies with an unsuppressed amplitude. In 
this paper we calculate this signature in detail. 

To circumvent the proposed test the modulus would 
have to decay before t ~ 1CU 22 s ~ which corre- 
sponds to m<j, = (m 2 ^) 1 / 3 - 10 12 GeV. 

A particularly interesting situation arises for super- 
string models with stabilized moduli that have at least 
one long lived modulus or modulus-like field whose mass 
is less than, or of the order of, the gravitino mass [BJ [7] . 
As mentioned above, the proposed test applies up to 
scalar masses of 10 12 GeV. Even though the tested mass 
range is well below the Planck scale, this would be a cru- 
cial test, because to our knowledge there is no other test 
for so massive, shortlived particles which are way beyond 
being testable in colliders, also in high precision tests. In 
this case the gravitino masses up to the same order of 
magnitude, 10 12 GeV, are probed. As far as we know, 
no other possibility has been proposed to probe so high 
supersymmetry breaking scales, albeit indirectly. 

Inflation does not only solve the horizon and flatness 
problems, but it also generates a scale invariant spectrum 
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of scalar and tensor (gravitational-waves) fluctuations, 
see e.g. [8]. Even though we know that the Universe was 
radiation dominated at BBN, its evolution history before 
that is unknown. 

In the following we examine the impact of a moduli 
dominated phase on the inflationary gravitational-wave 
background and discuss the prospects for observations in 
the inflationary gravitational-wave background. 



II. EVOLUTION OF GRAVITATIONAL WAVES 
FROM INFLATION 



Gravitational waves are the tensor perturbations h%j of 
the space-time metric, 

ds 2 = a 2 {7-j){-dif + (Sij + 2h ij )dx i dx j ) , (I) 

where a(j]) denotes the cosmic scale factor. The per- 
turbation hij is traceless, h\ — 0, and divergence free, 
d l hij = 0. The conformal time r\ is defined by dr\ — 
dt/a(t), where t denotes the physical time. The energy 
density of gravitational waves is then given by [HI [TP] 
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where G is Newton's constant. In this work an overdot in- 
dicates the derivative w.r.t conformal time 77. In Fourier 
space, the evolution of a gravitational-wave mode h in 
a Friedmann universe (neglecting anisotropic stresses) is 
determined by [5] 



h + 2-h + k 2 h = 0. 



(3) 



Introducing x = kr], and assuming power law expansion, 
a oc 77', this equation has the simple general solution 
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where j n and y n denote the spherical Bessel functions 
of order n as defined, e.g., in One might replace 

x/a by x x ~ q and adjust the pre-factors correspondingly. 
With this replacement it becomes evident that on super- 
Hubble scales, x < 1, the j-mode is constant while the 
y-mode behaves as x~ 2q+1 . From this general solution 
together with Eq. |2| one infers that p gw oc a~ 4 as soon as 
the wavelength is sub-Hubble, x > 1. (On super Hubble 
scales the 'energy density' of a mode is not a meaningful 
concept.) 

It is reasonable to assume that both modes have similar 
amplitudes after inflation, where x -C 1 for all modes of 
interest. If q > 1/2, the y-mode is decaying and soon 
after inflation we may approximate the solution by the 
j-mode. Note that a constant value of q corresponds to a 
constant background equation of state with the ratio of 
pressure P to energy density p given by 

w = P/p and q = 2/(3w + 1) . 



For a non- inflating (3w + 1 > 0) universe, q > 1/2 cor- 
responds to w < 1 and comprises all cases of interest. 
During inflation —1/3 > w > — 1 and q < — 1. 

In standard cosmology, the Universe is radiation dom- 
inated after reheating until the time of equality and mat- 
ter dominated afterwards. Therefore q = 1 until equality, 
where p ra( j = p ma ,t, and q = 2 after that. Since the en- 
ergy density in gravitational waves scales like radiation, 
its fraction is constant on scales which enter the hori- 
zon during the radiation dominated era and scales like 
a(rjk) oc rj 2 cx 1/k 2 for scales which enter during the mat- 
ter dominated era. A good approximation to the transfer 
function T 2 q (k) which relates the energy density per log- 
arithmic A:-interval to the amplitude of the gravitational- 
wave spectrum after inflation in standard cosmology is 
given in [12] . With this we obtain (for simplicity we ne- 
glect changes in the number of effective degrees of free- 
dom and the minor effect of today's vacuum domination) 



ri gw (fc) 
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?w(fc) =a.ad^T 2 q (fc), (5) 



where 



TcJk) = (1 + 1.57r; eq fc + 3.42( ?7eq fc) 2 )(r 7cmb fc)- 
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Here A-n is the amplitude of density fluctuations from 
inflation as measured in the CMB by the WMAP exper- 
iment P]. A|~2x 1CT 9 , and f2 rad ~ 5 x 1CT 5 . The 
ratio r is the tensor to scalar ratio which depends on the 
inflationary model, 77 oq and r? cm b are the conformal time 
at matter-radiation equality and at CMB decoupling re- 
spectively. This standard spectrum is indicated by the 
dotted line in Fig. [l] For Eq. ([6| a Harrison-Zel'dovich 
spectrum is assumed. For different primordial spectra 
with spectral index n s =/= 1 and ut 7^ of the primordial 
scalar and tensor fluctuations from inflation, the result 
^ has to be multiplied by (k/k c ) nT if the amplitude An 
and r (which then are scale dependent) are determined 
at the pivot scale k c . Changes in the background by non- 
standard evolution of the Universe have previously been 
studied in [14| . 

We now consider the situation that is likely to occur 
in superstring models, where soon after inflation moduli 
come to dominate when p^ — p rll d- We denote the cor- 
responding (conformal) time by i]^. We assume that the 
moduli decay briefly before nucleosynthesis at time t] c . 
We then compute the final gravitational-wave spectrum 
by matching the radiation solution (q = 1) before r?b to 
the matter solution (q = 2) at tj^ and back to the radi- 
ation solution at rj e . After ry e the Universe follows the 
standard evolution, so that the resulting spectrum sim- 
ply has to be multiplied by the standard transfer func- 
tion T 2 q (k). The generic shape of the resulting transfer 
function T is clear from the general solution: On super- 
Hubble scales the solution remains constant and T = 1. 
Scales that enter the horizon during the matter domi- 
nated phase at r)t> < % = 1/k < rj e are suppressed by 
a factor a(r)k) / a(i] e ) = (krj e )^ 2 since pgw & while 
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Pmat oc a 3 . Scales which have already entered before 
matter domination are maximally suppressed by a factor 
a(Vb)/a{r] c ) = (W^c) 2 - 

For sufficiently long matter domination, rj c /rj\, > 4, we 
find the following simple and accurate analytic approxi- 
mation to the exact result for the transfer function of an 
intermediate matter dominated phase: 



The gravitational wave spectrum today 
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where the best-fit analysis gives c = 0.5. The currently 
observable gravitational- wave spectrum is then simply 



, ad ^T e 2 q (fc)T 2 (fc; 7^,776) 



(8) 



with the fitting formula for T(k; rj e , r/b) from Eq. Q. The 
time when matter domination begins, r/b, is determined 
by the mass of the modulus field. Up to a coupling con- 
stant of order one, the mass also determines the time 
of the modulus decay, ?y c , which is the time when mat- 
ter domination ends. Using the general formula |15j for 
conformal time, 



T] = 1.5 x 10 5 
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-and assuming a large enough reheating temperature 
7t> ~ \J flbWplj ~ TtiQ- we find for a moduli mass 
of rriA ~ 100 TeV a value of % ~ 10~ 5 s . The resulting 
gravitational- wave background for rj c ~ ?7bbn as expected 
from the usual moduli problem is indicated as the thick 
solid line in Fig. [I] 

We have also indicated in Fig. [T] the results for moduli 
decaying much earlier, namely at T£ ~ 10 6 GeV (dashed 
line). Furthermore, for illustration we have indicated the 
suppression of the gravitational-wave spectrum in stan- 
dard cosmology from a particle with 30 TeV mass that 
enters thermal equilibrium after inflation and decays be- 
fore weakly-interacting massive particles (WIMP) freeze- 
out (thin solid line). In standard cosmology such signals 
might not only arise from any extension of the standard 
model containing long-lived particles like the axino or 
modulinos, but could even be desired [16]. 

We compare the spectra with present [T7] and fu- 
ture [18] gravitational-wave experiments in Fig. [5J see 
footnotes and references therein for details on the exper- 
iments. It is important that the decay temperature even 
of a much earlier modulus decay might be inferred from 
observation. To compare to the signal of an intermedi- 
ate reheating temperature see (19) . For larger moduli 
masses the moduli dynamics might depend directly on 
the reheating temperature. For the simplest potential, 
V = rn|0 2 /2, we estimate Tr ~ 10 9 GeV as the minimal 
reheating temperature corresponding to the sensitivity of 
BBO for our test to apply. 

Interestingly, pulsar timing arrays probe the scale of 
BBN, where the Universe is surely radiation dominated. 
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Figure 1: Intensity of inflationary gravitational waves vs. 
their frequency / = k/(2n) observed today. The thick solid 
line shows the expectation from the usual cosmological mod- 
uli problem. The dashed line corresponds to the case of a 
much earlier modulus decay at = 10 6 GeV. For compar- 
ison the dotted line shows a perfectly flat spectrum. The 
thin solid line demonstrates the impact on the spectrum of 
a particle with 30 TeV mass that enters thermal equilibrium 
after inflation and decays before WIMP freeze-out. Various 
frequencies corresponding to important and suggestive scales 
are highlighted by the vertical dashed lines: CMB indicates 
the scale of best sensitivity of CMB experiments, and the 
other frequencies relate to the horizon scale at the indicated 
event. 



Observation opportunities 
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Figure 2: Observation opportunities for the spectra of Fig.[T] 
In addition to the spectra of Fig. [T] sensitivity curves of exist- 
ing [T7j and future [18] gravitational-wave observatories are 
plotted. 



It would be particularly interesting, if they became sen- 
sitive to inflationary gravitational waves. The present 
CMB limit on r from gU] [H] is r < 0.2 on CMB scales. 
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III. CONCLUSIONS AND OUTLOOK 

We have shown that the gravitational-wave back- 
ground from inflation is strongly modified by an inter- 
mediate matter dominated era. If string compactifica- 
tion leads to an intermediate matter dominated phase, 
this will be observable in the gravitational-wave back- 
ground: on CMB scales one will detect the unmodified 
background from inflation, e.g., with the Planck satel- 
lite [22] or with a future CMB polarimeter [53]. However, 
on higher frequencies like milli-Hz or Hz probed by the 
gravitational- wave detectors indicated in Fig. [2] no signal 
will be detected. In other words, if these experiments 
will detect the signal from the inflationary gravitational- 
wave background as expected from the CMB, this will 
rule out all string compactifications that contain at least 
one scalar with a mass < 10 12 GeV — corresponding to 
the sensitivity of BBO — that acquires a large initial os- 
cillation amplitude after inflation and has only gravita- 
tional interaction strength. A correspondingly high su- 
persymmetry breaking scale, for example, of the order of 
the GUT scale may well render superstring models un- 
observable. 

Even though our derivation presented here is for an in- 
termediate matter dominated phase, the qualitative re- 
sult remains true also for a phase of thermal inflation. 
Such a phase would dilute gravitational waves on sub- 
Hubble scales even more strongly, and would render them 
undetectable for the experiments indicated in Fig. [2] 
while not affecting CMB scales. A (possibly) observable 
gravitational-wave spectrum created after thermal infla- 
tion were easily distinguishable from the primordial one 
by its shape [I]. 

Furthermore, from our derivation it is clear that for 
some other, non-inflationary, intermediate epoch starting 
at time r)b and ending at r] e , with equation of state P = 
wp with —1/3 < w < 1, the inflationary gravitational- 
wave spectrum will be suppressed (or enhanced for w > 
1/3) by a factor a(k) with 

n h g l al = a(k) x ttg with 

( 1 if k < Va 1 

a(k) (fc^^-D/^+D if v ^<k< v -\ 10) 

) / \2(3™-l)/(3w+l) K ' 



Our fitting formula for T(k) reproduces this behavior for 
an intermediate matter dominated era, w = 0. For ic/0 
the exponents ±2 in the denominator would have to be 
replaced by ±2(1 - 3iu)/(l + 3iu). 



Of course there is the caveat that (as also many 
other inflationary models) inflation in models inspired 
by string theory [21], depending on the compactification 
and moduli stabilization, may predict only a very low 
gravitational-wave background that cannot be measured 
by proposed experiments, neither in the CMB nor di- 
rectly on smaller scales. In this case, an experiment that 
would be able to detect the background generated in the 
intermediate matter dominated era as predicted in |25j 
should be conceived. Such a background is, however, sup- 
pressed with respect to the amplitude from inflation by 
the ratio of the corresponding Hubble rates (ffb/#inf) 2 , 
where denotes the Hubble rate at the beginning of 
the intermediate matter dominated phase and -ffinf the 
one during inflation. 



In conclusion, we have shown that combining fu- 
ture CMB polarization measurements with very sensi- 
tive gravitational-wave probes such as BBO can provide 
a crucial test for a large class of string models. Since 
string theory is our best candidate for a fundamental 
theory of Nature, and there has not been proposed any 
experiment to test it, it is of uttermost importance that 
we realize a BBO-like experiment. 
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